Persistent holes in a fluid 



m 

O 

o 

(N 

> 
O 



> 

O 



o 



> 
X 



Florian S. Merkt, Robert D. Deegan, Daniel I. Goldman, Erin C. Rericha, and Harry L. Swinney 
Center for Nonlinear Dynamics, The University of Texas at Austin, Austin, TX 78712 

(Dated: February 8, 2008) 

We observe stable holes in a vertically oscillated 0.5 cm deep aqueous suspension of cornstarch 
for accelerations a above lOg. Holes appear only if a finite perturbation is applied to the layer. 
Holes are circular and approximately 0.5 cm wide, and can persist for more than 10^ cycles. Above 
a ~ llg the rim of the hole becomes unstable producing finger-like protrusions or hole division. At 
higher acceleration, the hole delocalizes, growing to cover the entire surface with erratic undulations. 
We find similar behavior in an aqueous suspension of glass microspheres. 



The free surface of a fluid at rest in a container is flat. 
Departures from flatness induce a restoring flow whether 
the fluid is Newtonian, viscoelastic, or liquid crystalline: 
poke the surface, and the resulting indentation will be 
filled by the ensuing flow. In contrast, we have discovered 
that a vibrated aqueous suspension of cornstarch or glass 
microspheres can permanently support holes and vertical 
finger- like protrusions. 

The catalog of interface morphologies in accelerated 
fluids is broad and well documented. Sinusoidal accel- 
eration produces Faraday waves dj, solitons and 
jets Impulsive acceleration produces the Richtmyer- 
Meshkov instability |^, and continuous acceleration 
the Rayleigh- Taylor l7| instability, which are character- 
ized by spires and bubbles The holes and fingers 
in vibrated cornstarch are distinct from these phenom- 
ena because they display neither the oscillation about 
a flat state of Faraday waves and solitons, the finite life- 
time of jets, nor the unbounded growth of the Richtmyer- 
Meshkov or Rayleigh- Taylor instabilities. 

Figure n shows holes in a vibrated aqueous suspension 
of cornstarch (a-b) and glass microspheres (d-e). De- 
pending on the container acceleration and frequency, the 
initial fluid surface is either flat or corrugated by Fara- 
day waves. Above a critical acceleration, a finite localized 
perturbation of the fluid surface will grow into a stable 
cylindrical void that extends from the top to the bottom 
of the fluid layer. The most noteworthy feature of holes 
is their permanence: they do not close despite the hy- 
drostatic pressure of the surrounding fluid, persisting as 
long as our observation (> 10^ cycles). At yet higher 
accelerations, the holes lose stability to a temporally and 
spatially erratic state (see Fig.^c) & (f)). 

Experiment. A layer of cornstarch or glass micro- 
spheres in liquid was vertically oscillated sinusoidally 
with a frequency / from 50 to 180 Hz and a peak ac- 
celeration a up to 275. The fluid container was round 
with a 9.4 cm diameter, and made from an aluminium 
base plate and acrylic sidewalk The container was at- 
tached to the shaker through an insulating rod to avoid 
heat transfer from the shaker to the container, and was 
sealed with an acrylic top to reduce evaporation. The 
layer depth was 0.5 cm for the cornstarch mixture and 
0.2 cm for the glass microspheres. The peak acceleration 



of the container was controlled to within ztO.Olg. 

The patterns were recorded with a 30 frames per sec- 
ond digital camera and a 1000 frame per second video 
camera. Lighting was provided by a ring of LEDs strobed 
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FIG. 1: Top view of a vibrated layer of (a-c) cornstarch 
and of (d-f) glass microspheres in liquids matched in density 
to the particles. Diameter is 9.4 cm. White corresponds to 
the highest points, and black to depressions that reach the 
container bottom. Holes without Faraday waves: (a) a = 12g, 
f = 150 Hz; (d) 30g, 100 Hz. Holes with Faraday waves: (b) 
123, 60 Hz; (e) 27.83, 92 Hz. Delocalized state: (c) 30g, 
120 Hz; (f) 30g, 60 Hz. 
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at / or //2. The illumination was such that peaks on the 
surface appear bright and valleys dark. 

We used unmodified, regular cornstarch consisting of 
27% amylose and 73% amylopectin from Sigma Aldrich. 
The powder was dried at 50° C for a week and stored in a 
desiccator. A mixture was prepared daily by combining 
23.22 g of cornstarch, and 36.78 g of a density-matched 
aqueous solution of CsCl with a density of 1.68 g/cm^. 
The viscosity of the suspension was measured in a plate- 
plate geometry with a Paar Physica TEK150 rheometer. 
Though all the quantitative results presented here are for 
cornstarch, we also used a glass microsphere and water 
mixture for the qualitative comparison in Fig. ^ The 
glass mixture was prepared with 29.37 g of 1 — 20 /im 
glass balls from Jaygo (Union, NJ), and 18.36 g of an 
aqueous solution of sodium polytungstate with a density 
of 2.55 g/cm^. 

Hole dynamics and size selection. The interface of a 
vibrated fluid is flat at low accelerations and rippled by 
Faraday waves above a critical acceleration [jj]. Holes 
form in the cornstarch mixture only if an indentation 
deeper than about 50% the layer depth and 0.4-2.0 cm 
wide is applied to the surface. We used a puff of com- 
pressed air delivered through a nozzle pointed at the sur- 
face to form a crater. The subsequent evolution of the 
crater depends on the parameters a and /. For all param- 
eter values explored, the crater may close within seconds. 
In the region of the phase diagram (Figl^J marked "un- 
stable" all craters are annihilated by this rapid decay. 
In regions marked "stable" and "metastable" , the crater 
may develop into a hole: the crater adjusts to a preferred 
size (as shown in FigOJ and remains fixed at that size for 
an extended period of time. At / = 120 Hz the charac- 
teristic diameter of the holes is 0.4 cm and at / = 180 Hz 
0.6 cm. Holes have a broad distribution of lifetimes, even 
for fixed / and a; in the metastable phase the maximum 
lifetime is less than 10^ cycles and in the stable phase it 
is greater than 10^ cycles. 

Holes collapse via one of two mechanism depending on 
their lifetime. Short lived holes collapse within a few 
seconds by a uniform contraction. Longer lived holes 
with lifetimes of order lO'* holes develop a hump on their 
rims, which then falls onto the hole either covering it or 
rendering it so small that it collapses rapidly by uniform 
contraction. On occasion this latter mechanism causes 
the hole to divide rather than collapse. 

The cross-section of a hole is typically circular with a 
slight eccentricity. The shaft is smooth and slightly ta- 
pered towards the base. At the base of the hole, there is 
usually a ribbon of material about 0.1 cm high and wide 
that bisects the shaft. These features are not visible in 
Fig. m due to the lighting. At accelerations immediately 
above the stability line the material surrounding a hole 
remains level with the rest of the layer. At higher ac- 
celeration, the hole's rim rises about 0.05 cm above the 
surrounding fluid. 
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FIG. 2: Phase diagram for a vibrated aqueous cornstarch 
suspension as a function of acceleration and frequency. For 
accelerations above the o, a surface perturbation will form a 
hole. In the metastable regime holes always collapse within 
10^ container oscillations; in the stable regime some holes 
endure for at least 10^ oscillations, show the acceleration 
threshold, below which holes collapse in less than 10* cycles. 
In the delocalized regime the state is chaotic as shown in 
Fig.n^c). Faraday waves appear with increasing acceleration 
at and disappear with decreasing acceleration at T. 



In addition to a slow meandering over hundreds of 
oscillations, hole size oscillates synchronously with the 
driving frequency, as shown in the inset of Fig. |3 The 
radius of a hole evolves as r{t) = r + Sr sin (27r/f -I- 9), 
where r is the mean radius, Sr is the amplitude of the 
oscillation, and 9 is the phase relative to the container 
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FIG. 3: Hole radius versus time in the cornstarch suspension 
for 6 different trials at 15g and 150 Hz. At t = a jet of 
gas was applied to the surface for 2 s. The hole grows during 
this forcing, and immediately begins to shrink after the gas is 
cutoff. Inset: rapid oscillation of the hole. Note time scale. 



motion z{t) = sin (27r/t) where Zo = — a/(27r/)^ and 
z is positive when the container is above its rest point. 
Sr is typically 0.15r and also varies by about 20% on a 
timcscale of order 10'^/"^ . The phase lag is shown as a 
function of frequency in Fig. 0] 
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FIG. 4: The tangent of the phase lag versus frequency at 
a = 15g The solid line is a fit to — {f / foY) where a and 
/o are fitting parameters. 



Holes do not interact when their centers are separated 
by more than two diameters. Therefore, as shown in 
Figure^a), holes do not form regular patterns. Further- 
more, they can be located anywhere in the container. 
Holes do occasionally come sufficiently close to interact, 
and do so by merging or one causing the other to collapse. 

Faraday waves. Though Faraday waves are excited in 
cornstarch suspensions, the formation of holes is unre- 
lated to these waves, as can be seen from the existence of 
holes when no Faraday waves are present, as in Fig. 1(a) 
and (d). Furthermore, the phase boundaries for holes 
and Faraday waves are distinct. The Faraday transition 
is hysteretic and, moreover, the flat and surface wave 
state can coexist. As acceleration is increased at a fixed 
frequency, small patches of surface waves appear. The 
onset of these are indicated by the A in Fig. |5] As the 
acceleration is raised further, the patches grow and ul- 
timately engulf the entire surface at values typical 50% 
higher than onset. Upon decreasing the acceleration, the 
surface waves breakup into patches which finally disap- 
pear at accelerations indicated by the T in Fig.|2| Despite 
the ambiguity, there is a clear difference in the param- 
eter dependence for the onset of holes and the onset of 
Faraday waves. 

Delocalization. In the "delocalized" region in Fig. 12 
a perturbation will generate a hole which immediately 
grows a protrusion from its rim. The sequence of photos 
in Fig. El shows the growth of such a structure. They are 
surprisingly tall and long-lived, reaching heights of 2 cm 
and remaining upright for thousands of oscillations. The 
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FIG. 5: Side view of the first steps of a delocalizing hole 
in cornstarch. These photographs were taken every 0.9 s. 
Time increases from left to right and top to bottom. An 
initial hump on the rim begins growing upward, it reaches 
a maximum height, and then topples outward enlarging the 
area of fluid motion. This process is repeated until the entire 
surface of the liquid is active in the creation and destruction 
of vertical structures and voids (see Fig. 0^c)). These images 
were taken at a —25g and f — 80 Hz. 



protuberance eventually falls, nucleating a new hole, and 
the growth of a finger is repeated. After many repeti- 
tions the entire surface writhes with fingers and holes 
yielding the spatially and temporally erratic state shown 
in Fig.mc) and (f). 

Shear Thickening. We attribute the formation of holes 
to shear thickening, i.e., an increase of viscosity with 
shear rate. A plot of the viscosity for two identically 
prepared cornstarch suspensions is shown in Fig. El As is 
typical for shear thickening fluid 10], an initial decrease 
of the viscosity for increasing shear rate is followed by a 
rapid increase at a critical shear rate, which in this case 
is 7c = 8 s~^. Our conclusion is predicated on four obser- 
vations: the critical shear rate is similar in magnitude to 
the shear rate at the throat of a hole; the phase lag of the 
holes radius relative to the driving force indicates the ma- 
terial response is primarily viscous; holes are unstable in 
a non-shear thickening fluid; and holes are stable in shear 
thickening fluids other than cornstarch suspensions. 

The size of the hole varies synchronously with the 
driving frequency. The root-mean-squared shear rate at 
throat of the hole is approximately the interface speed, 
(27r/\/2)(5r/, divided by a length scale of order the depth 
of the layer h. Using Sr/r ~ 0.15, / = 100 Hz, and 
r/h ~ 0.4 yields a shear rate 7 ~ 27 s^^. The similarity 
of this value to the critical shear rate strongly suggests 
that shear thickening is an essential ingredient for holes. 

The frequency dependence of the phase lag 6 indi- 
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FIG. 6: Apparent viscosity versus maximum shear rate ')max 
for two identical samples of an aqueous cornstarch suspension. 
The torque r on the plate was measured using a plate-plate 
geometry with a gap spacing h of 0.10 cm and a plate radius 
R of 2.50 cm. The apparent viscosity was calculated as it 
would for a Newtonian fluid: r) = 2t / {-K'^rnaxR^) ■ Since the 
shear rate increases with the radius, our measurement of the 
viscosity is a convolution of the geometry and the response of 
the fluid. The dotted line represents measured values that are 
not reliable due to an instability in the flow. The cornstarch 
samples were the same as in Fig. and the temperature was 
kept at 27.6 ± 0.2 °C. 



cates a combined viscous and elastic response of the fluid. 
Modelling the hole motion as a spring in parallel with a 
dashpot (i.e., a Voigt element, see for example) yields 
tan^ = a/{l — {f/fo)^), where a is a constant propor- 
tional to the dissipation and fo is the resonant frequency. 
As shown in Fig. ^ the data are well modelled by this 
equation;the resonance frequency is around 58 Hz. Since 
the Voigt element is primarily viscous at high frequency, 
the phase data indicates that at frequencies above 60 Hz 
the material response becomes dominantly viscous; it is 
noteworthy that stable holes only form in this high fre- 
quency region. 

We also tried to form holes in silicon oil (a Newtonian 
fluid) and various polybutadienes with viscosities 8.3, 28, 
and 72 Poise. In all cases, however, the initial crater was 
backfilled within a few hundred container oscillations. 

In Barnes' review of shear thickening fluids, he asserts 
that any sufficiently concentrated suspension of solids in 
a fluid will shear thicken If our hypothesis that 



shear thickening is crucial for hole formation is correct, 
then any sufficiently dense suspension of solids ought to 
support persistent holes. For a dense suspension of glass 
microspheres this is indeed the shown in Fig.^d) 

&(e). 

Conclusions. We have shown that a vertically vi- 
brated aqueous suspension of cornstarch displays a num- 
ber of unexpected patterns-holes, fingers, and a delocal- 
ized state-that can only be attained by the application of 
a finite perturbation. Holes have a preferred size, vibrate 
synchronously with the container, and are unrelated to 
the well-known Faraday wave instability. At higher ac- 
celeration the hole becomes unstable to the creation of 
finger-like structures, which at even higher accelerations 
are responsible for the transition to a delocalized state in 
which the free surface is covered with fingers and holes. 
Our observation suggest that shear thickening is the key 
ingredient for hole formation. The outstanding questions 
are how do holes stay open, how can the finger like pro- 
trusions grow and remain unright, and what drives the 
transition to the delocalized state? 

We thank M. Schroter, W.D. McCormick, P. HeU, 
A. Lee, W. Zhang for valuable suggestions, and P. Green 
and A. Jamieson for assistance with the viscosity mea- 
surements. We acknowledge support by the Engineering 
Research Program of the Office of Basic Energy Sciences 
of the U.S. Department of Energy under grant number 
DE-FG03-93ER14312. FSM acknowledges support by 
the Fricdrich Ebert Stiftung. 



[1] M. Faraday, Philos. trans. R. Soc. London 121, 319 
(1831). 

[2] O. Lioubashevski, H. Arbell, and J. Fineberg, Phys. Rev. 

Lett. 76 3959 (1996). 
[3] O. Lioubashevski, Y. Hamiel, A. Agnon, Z. Reches, and 

J. Fineberg, Phys. Rev. Lett. 83, 3190 (1999). 
[4] M.S. Longuet-Higgins, J. Fluid Mech.127, 103 (1983). 
[5] R.D. Richtmyer, Commun. Pure Appl. Math. 13, 297 

(1960). 

[6] E.E. Meshkov, Fluid Dyn. 4, 101 (1969). 

[7] G.L Taylor, Proc. Royal Society, A 201, 192 (1950). 

[8] U. Alon, J. Hecht, D. Ofer, and D. Shvarts, Phys. Rev. 

Lett. 74, 534 (1995). 
[9] J. D. Ferry, Viscoelastic Properties of Polymers (Wiley, 

1961). 

[10] H.A. Barnes, J. Rheol. 33, 329 (1989). 



